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1 Introduction 



The definition of fundamental objects appearing in perturbative QCD calculations, such 
as the color charge and parton distribution functions, depends on some unphysical quan- 
tities. Theoretical predictions for physical quantities are independent of these unphysical 
quantities if the appropriate power expansions are summed to all orders. However, in 
practice, we are able to perform only finite order calculations, which provide theoretical 
predictions depending on the choice of the numerical values of these unphysical quanti- 
ties. The convenient choice of these numerical values thus plays very important role in 
finite order calculations. There are three general methods how to fix the renormalization 
scale: the Principle of Minimal Sensitivity (PMS) [1], the Effective Charges (EC) [2] and 
Brodsky-Lepage-MacKenzie (BLM) approach [3]. They all recognize the fact that the ex- 
istence of well defined natural physical scale of a given hard process (like in DIS) does 
not by itself imply that the renormalization scale should be identified with it, but start 
from different strategies how best to fix it. PMS looks for the region of local stability of 
finite order perturbation calculation with respect to the variation of the renormalization 
scale, EC (sometimes also called Fastest Apparent Convergence, FAC approach) prefers 
the value where all higher order contributions vanish, and BLM follows closely the recipe 
used in QED. The PMS and EC approaches can be applied also for fixing the factoriza- 
tion scale. All three approaches, together with the conventional approach in which the 
renormalization scale is identified with the natural scale characterizing the hardness of the 
process and allowed to vary within some reasonable range [4], have been extensively used 
in phenomenological analyses of hard scattering processes. See, for example, ref. [5] for 
extensive and detailed comparison of these approaches in QCD analysis of event shapes 
measured at LEP. There is also alternative approach to the formulation of perturbative 
QCD, which relates directly predictions for different physical quantities [6, 7]. However, 
little attention has so far been paid to the freedom in the choice of the factorization scheme 
even though it is, in principle, as important as the choice of the factorization scale. All 
phenomenological analyses of hard scattering processes carried out so far have been per- 
formed in the MS and DIS factorization scheme. This is probably related to the fact that 
the freedom in the choice of the factorization scheme is enormous, even at the NLO. 

The aim of this paper is the investigation of the freedom associated with the definition 
of parton distribution functions. The immediate motivation for this study is the potential 
exploitation of the freedom in the choice of the factorization scheme for the construction 
of consistent NLO Monte Carlo event generators [8-10]. 

The paper is organized as follows. The next section contains a review of basic facts, 
relations and notation. Section 3 is devoted to the discussion of the freedom associated 
with the factorization procedure in massless perturbative QCD. The subject of this section 
is the characterization of the freedom associated with factorization, the specification of 
factorization schemes via the corresponding higher order splitting functions and an overview 
of general formulae which govern the dependence of parton distribution functions and 
hard scattering cross-sections on unphysical quantities associated with the factorization 
procedure. Their application at the NLO is described in Section 4. Since the so called 
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ZERO factorization scheme appears as the optimal factorization scheme for NLO Monte 
Carlo event generators, Section 5 is devoted to the numerical analysis of its practical 
applicability at the NLO. The summary and conclusion are presented in Section 6. Some 
important technical details are the subject of appendices. 

2 Basic facts and notation 

Within the framework of perturbative QCD, theoretical predictions for physical quantities 
are calculated as power expansions in the renormalized coupling parameter of QCD a(/x, RS) 
(a = Us/ it). In the case of massless perturbative QCD, which will be the object of our 
interest, the renormalized coupling parameter a(/x, RS) is the only free parameter that 
characterizes the theory.^ However, the renormalized coupling parameter a(/i, RS) is not 
a single number, but it is a function of the so called renormalization scale ^ and the 
parameters that specify the renormalization scheme RS. Both the renormalization scale 
and the renormalization scheme fix the ambiguities associated with the renormalization 
procedure, which removes ultra violet singularities from perturbative calculations.^ The 
dependence of the renormalized coupling parameter a(;U, RS) on the renormalization scale 
// is determined by the differential equation 



where b = (33 — 2nf)/6, c = (153 — 19nf)/(66 — 4nf) and rif denotes the number of 
quark flavours. Whereas the first two coefficients are unique, the higher order coefficients 
q(RS) are completely arbitrary numbers and, together with the initial condition of the 
preceding differential equation, can be used for the unique specification of the corresponding 
renormalization scheme RS. 

The coefficients of the power expansions that represent theoretical predictions for phys- 
ical quantities depend on the renormalization scale and the renormalization scheme in such 
a way that if the power expansions are summed to all orders, then the obtained theoretical 
predictions are independent of the renormalization scale and the renormalization scheme. 
However, the finite order theoretical predictions which we obtain by truncating the corre- 
sponding power expansions depend on these unphysical quantities. The convenient choice 
of the renormalization scale and the renormalization scheme thus plays very important role 
in finite order calculations^ (in practice, we are reliant only on finite order calculations). 

^If the number of quark flavours is given. 

^The ultra violet singularities are absorbed into the definition of the renormalized coupling parameter 
a(/i,RS). The ambiguity of the renormalization procedure arises from the fact that the renormalized 
coupling parameter a(/i, RS) is not defined uniquely by the requirement of the absorption of the ultra violet 
singularities — an arbitrary finite term can be absorbed along with every singularity. 

^The ambiguity of the renormalization procedure does not necessarily have to be seen only as an incon- 
venience because a convenient choice of the renormalization scale and the renormalization scheme (which 
are not fixed but may depend on the process and its kinematics) can improve the agreement between finite 
order theoretical predictions and experiment. 



da(|U, RS) 



dln^ 
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To obtain theoretical predictions for processes whose initial state involves hadrons, 
we need to know the parton distribution functions Di/f{{x), which describe the parton 
structure of the relevant hadrons.^ As a simple example, consider deep inelastic lepton- 
hadron scattering. In this case, the relevant cross-sections can be expressed in terms of 
structure functions F(^x,Q'^^, which, according to the factorization theorem, are given [11] 
as the convolution integral^ 

F{x,Q^)=Y, T — Qf-,Q',M,FS,RS) A(y,M,FS,RS) (2.2) 

^ Jx y \y J 

where (x, , M, FS, RS) stands for the corresponding coefficient functions and Di{x, M, FS, RS) 
represents the relevant parton distribution functions. Both the coefficient functions Cj (x, , M, FS, RS) 
and the parton distribution functions Di{x, M,FS,RS) depend on the factorization scale 
M, the factorization scheme FS and the renormalization scheme RS, which is used for 
the factorization procedure, but the structure function F(x,(5^) is independent of these 
unphysical quantities (at least if all relevant expansions are summed to all orders), which 
fix the ambiguity associated with the factorization procedure, which removes the so called 
collinear singularities from expressions for physical quantities.^ 

The parton distribution functions cannot be calculated from perturbative QCD and 
therefore must be taken from experimental data. Perturbative QCD determines only their 
dependence on unphysical quantities — the factorization scale, the factorization scheme 
and the renormalization scheme used for the factorization procedure. The dependence on 
the factorization scale is described by the evolution equations 

'''^^IffiM''''^ = ^ 7 (i' ^' ^'-'^ 

where the splitting functions Pjj(x, M, FS, RS) can be expanded in powers of a(M, RS) 

oo 

Pij{x, M,FS,RS) = ^a'=(M,RS)i^5;''^(x,FS,RS). (2.4) 

k=0 



Di/ h{x)Ax is the number of partons of species i inside hadron H each of which carries the fraction of 
the hadron momentum that is between x and x + Ax. If it is not necessary to specify to which hadron parton 
distribution functions belong, then the corresponding designation of the parton distribution functions will 
be dropped in the following text. 

^This formula represents a separation of short distance properties of the theory described by the coef- 
ficient functions Ci (x, Q^, M, FS, RS) and large distance properties of the theory described by the parton 
distribution functions Di(x, M, FS, RS). The validity of this formula is not constrained only to perturbative 
QCD. 

^The treatment of the factorization procedure in this text is based on perturbative calculations and 
therefore the renormalized coupling parameter that is used for these calculations has to be specified (the 
factorization procedure thus has to be preceded by the renormalization procedure). The used renormalized 
coupling parameter is specified by the factorization scale AI and the renormalization scheme RS. Within 
the framework of the factorization procedure, the collinear singularities are absorbed into the definition of 
the parton distribution functions. Analogously to the case of the definition of the renormalized coupling 
parameter, the (dressed/renormalized) parton distribution functions are not determined uniquely by the 
requirement of the absorption of the collinear singularities. The associated ambiguity is then fixed by the 
specification of the factorization scheme FS. 
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Whereas the LO sphtting functions P.^ (x) are unique (they are independent of the factor- 
ization scheme FS and the renormahzation scheme RS), the higher order sphtting functions 
pj^j\x, FS, RS), k > 1 are completely arbitrary functions and can be used for labelling fac- 
torization schemes. 

The coefficient functions Ci{x, , M, FS, RS) are, at least in principle, fully calculable 
within the framework of perturbative QCD and can thus be expanded in powers of an 
arbitrary renormalized coupling parameter: 

oo 

a {x, Q^, M, FS, RS) = a''(/^> ^Se) \x, Q^,fx, RSe, M, FS, RS) . (2.5) 

The coefficient functions Cj (x, Q^, M, FS, RS) are independent of the renormahzation scale 
/X and the renormahzation scheme RSe if the corresponding expansions are summed to all 
orders (but still depend on the factorization scale M, the factorization scheme FS and 
the renormahzation scheme RS). The renormahzation scheme RSe, which is employed for 
expanding the coefficient functions Cj (x, Q^, M, FS, RS) , is in principle different from the 
renormahzation scheme RS, which is used for the factorization procedure. However, in 
practice both renormahzation schemes are chosen to be identical, which simplifies calcu- 
lations. Further simplification of calculations can be achieved by setting fi = M, which 
is a common and legitimate choice, but there is a good reason for treating at least the 
renormahzation scale and the factorization scale M as independent of each other: the 
renormahzation scale emerges from the renormahzation procedure, which deals with ultra 
violet singularities, which are related to short distance properties of the theory, while the 
factorization scale appears in the factorization procedure, which treats collinear singulari- 
ties, which are connected with large distances. Keeping these scales separate can also lead 
to a better agreement between finite order theoretical predictions and experimental data. 

As a simple example illustrating the freedom in the choice of the factorization scheme, 
consider the structure function F-^Nsi.^^Q'^) which is defined as 

F,%{x,Q') = ^[F,^^x,Q')-Frix,Q')). (2.6) 
The structure function F2^g(^x,Q'^) can be approximately expressed as^ 

^2!ns(^'Q') =^'yC'Ns(^,Q',M,FS^9Ns(y,M,FS) (2.7) 
where (the arguments are suppressed for brevity) 

Cns = 3 (C2,u - C2,d), qm = 2 (-^^u/p + -Du/p - -Dd/p - -Dd/p) • (2-8) 

^In the following, we limit ourselves only to such factorization schemes whose definition preserves the 
symmetry between quarks and antiquarks and between quark flavours. For simplicity, the renormahzation 
scheme used for the factorization procedure is fixed and is the same as that in which the renormalized QCD 
coupling parameter is defined. Hence, the dependence on the renormahzation scheme is not written out 
explicitly. 
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The quark non-singlet distribution function q-t<{s{x, M,FS) satisfies tlie evolution equation 
d,.s(x,MFS) ^ r^P,sf^^,M,Fs),Ns(.,M,FS). (2.9) 

dinM Jx y \y J 

At the NLO approximation, we retain only the first two terms in the expansions of the 
coefficient and splitting function: 

CNs(a;,Q^M,FS) =5(l-x) + a(/i)ci,s^(2;,Q^M,FS), (2.10) 

Pns(x,M,FS) = P^l{x) + a{M)P^l{x,YS). (2.11) 

While i'Ns(x) is unique, both Pj^g'(x,FS) and Cj(jg^ (x, Q^, M, FS) are arbitrary, subjected 
only to the relation [12] 

cW(x,Q2,M,FS) =p(°)(x)ln^ + lp«(x,FS) + K(:E) (2.12) 

where k(x) is a factorization invariant, i.e. independent of both the factorization scale and 
the factorization scheme. This relation connects the coefficient function (x, Q'^,M, FS) 
associated with one-loop perturbative computations with the splitting function P^g{x, FS) 
associated with two-loop perturbative computations. This connection reflects the fact 
that the finite parts associated with one-loop perturbative computations, which cause the 
ambiguity of the coefficient function C^^, also appear in two- loop perturbative computation 
of P^^- However, relation (2.12) does not imply that the splitting function pj^^ can be 
determined from a one-loop perturbative calculation of the coefficient function C^^ because 
determining the factorization invariant k{x) requires a two-loop perturbative calculation. 
Relation (2.12) only expresses the fact that the ambiguity of the coefficient function 
is correlated with the ambiguity of the splitting function P^g- A more detailed discussion 
of relation (2.12) can be found, for instance, in [13]. 

The ambiguity associated with the factorization procedure is large, but almost unex- 
ploited in practice [12]. The most widely used factorization scheme is the so called MS 
factorization scheme, which is suitable for theoretical calculations. In this factorization 
scheme, both the splitting function P^^ (x, FS) and the coefficient function C^g^ (x, Q^, M, FS) 
are nonzero. Some analyses are also performed in the DIS factorization scheme, which is in- 
troduced in [14] in order to express the relation between the structure function F2^(^x, Q^) 
and parton distribution functions in the same way as in the parton model (this can be done 
for only one, but arbitrary, ratio of Q and M), which means that the coefficient function 
C^^g (x, Q^, M, FS) vanishes for M = Q. This implies that all NLO corrections are in- 
cluded in the NLO splitting function Pj^g''(x,FS) and thus exponentiated by the evolution 
equation. A certain opposite to the DIS factorization scheme is represented by the ZERO 
factorization scheme in which all NLO splitting functions vanish and therefore all NLO cor- 
rections are contained in the hard scattering cross-sections (in our case this means in the 
NLO coefficient function C^^(^x,Q'^,M,FS^), which implies that the NLO corrections are 
completely unexponentiated. The ZERO factorization scheme is mentioned, for instance. 
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in [13]. In the next paragraph, it is explained why the ZERO factorization scheme could 
be useful for constructing consistent NLO Monte Carlo event generators. 

At present time many QCD cross-sections at parton level are known at the NLO 
accuracy and necessary algorithms for their incorporation in Monte Carlo event generators 
have been developed. However, these algorithms attach to them initial state parton showers 
only at the LO accuracy because no satisfactory algorithm for generating initial state 
parton showers at the NLO accuracy (in the standard MS factorization scheme) has been 
found so far. Since the initial state parton showers induce the scale dependence of parton 
distribution functions, it is inconsistent to attach LO initial state parton showers to NLO 
QCD cross-sections, which include NLO parton distribution functions. This deficiency 
could be removed by exploiting the ZERO factorization scheme, in which the NLO initial 
state parton showers are formally identical to the LO ones. The main advantage of this 
approach is the fact that the existing algorithms for parton showering and for attaching 
parton showers to NLO cross-sections need not be changed. The only necessary action is to 
transform hard scattering cross-sections from the standard MS factorization scheme to the 
ZERO factorization scheme and to determine parton distribution functions in the ZERO 
factorization scheme. 

The construction of consistent NLO Monte Carlo event generators in which initial state 
parton showers can be taken formally at the LO thus constitutes one of the motivations 
for investigating factorization schemes. Hence, the following section will be devoted to a 
more detailed discussion of the freedom associated with factorization, with the emphasis 
on its quantification. 

3 The freedom associated with factorization in massless QCD 
3.1 Basic facts about factorization in massless QCD 

The parton model description of a hard hadronic collision has three basic ingredients: the 
parton distribution functions of the colliding hadrons, a set of parton cross-sections that de- 
scribes the hard scattering of the partons in the process and some model for the hadroniza- 
tion of the final state partons into observable hadrons (the description of hadronization is 
not necessary if the hadrons in the final state of the hard hadronic collision are not spec- 
ified). Parton cross-sections are, contrary to parton distribution functions and hadroniza- 
tion, fully calculable within the framework of perturbation theory. However, QCD radiative 
corrections to parton cross-sections are typically divergent (even after the renormalization 
procedure because the appropriate singularities are related to the behavior of massless 
perturbative QCD at large distances) and therefore useless for the straightforward incor- 
poration in the parton model. Fortunately, according to the factorization theorem [15], the 
singularities connected with the partons in the initial state are process independent and can 
be extracted from the parton cross-sections and absorbed into the bare parton distribution 
functions Diix) of the naive parton model. ^ The convolution of the singular factors from 

*The elimination of other possible singularities from parton cross-sections, which is related to the defi- 
nition of the final state, is independent of the factorization of the singularities connected with the partons 
in the initial state. 
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parton cross-sections with the bare parton distribution functions, denoted by Di{x) in the 
following, then defines the dressed parton distribution functions M, FS, RS), which 

are finite, physically measurable, process independent, but ambiguous because arbitrary 
finite terms can be added to the singular terms that are absorbed into the definition of the 
dressed parton distribution functions. The ambiguity associated with the definition of the 
dressed parton distribution functions M, FS, RS) is discussed below. 

Within the framework of dimensional regularization in d = 4 — 2e space-time dimen- 
sions, the relation between dressed and bare parton distribution functions is given by the 
convolution ^ 

A(x,M,FS,RS) = ^ / ^A,/-,M,FS,RS] (3.1) 

where the functions M, FS, RS), which represent the singular factors that are ab- 

sorbed into parton distribution functions, can be expanded in powers of a(M, RS): 

oo 

Aij{x, M,FS,RS) = ^a^(M,RS)^Sf (a;,FS,RS), FS, RS) = 5,j6{l - x) (3.2) 

and the higher order coefficients of the preceding expansion can be expressed in the form 

oo ^ 

^Sf (x,FS,RS) = ^^ylSf (2;,FS,RS), k>l (3.3) 

1=0 

where the functions A^^'''(x, FS, RS) are independent of e. The functions A^j^\x), which 
represent the finite part of the absorbed factors, can be chosen arbitrarily and their choice 
defines the factorization scheme FS. The functions ^^^^^(x, FS, RS), / > 1, which specify 
the singular part of the absorbed factors, are then uniquely determined by the properties 
of the theory and the choice of the factorization scheme FS (the choice of the functions 

A detailed analysis of the freedom associated with factorization in massless perturba- 
tive QCD is presented in Appendix D. The following two subsections contain an overview of 
important results of this analysis. The rest of this paper, with the exception of Appendices 
B and D, is concerned only with the case of four space-time dimensions {e = 0). Some 
relations are expressed in terms of Mellin moments, which are defined in Appendix A. The 
last subsection of this section is devoted to a general discussion of the issue of practical 
applicability of factorization schemes that are specified by the corresponding higher order 
splitting functions. 

3.2 Parton distribution functions 

In this subsection, we present formulae describing the change of the unphysical quantities 
on which parton distribution functions depend and some important facts concerning the 
freedom associated with the factorization procedure. 

Changing the factorization scheme in which parton distribution functions are defined 
is described by the formula 

D(n,M,FSi,RS) = T(n, M, FSi, FSa, RS)D(n, M, FS2, RS) (3.4) 
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where the multiphcation is matrix multiphcation. The parton distribution functions are 
represented by a column vector and the transformation matrix T(n, M, FSi, FS2, RS) is a 
square matrix. The preceding formula determines the change of the factorization scheme 
for the fixed factorization scale M and renormalization scheme RS, which is used for the fac- 
torization procedure. The transformation matrix T(n, M, FSi, FS2, RS) can be expanded 
in powers of a(M, RS): 

00 

T(n,M,FSi,FS2,RS) = ^ a'=(M, RS) T('^)(n, FSi, FS2) (3.5) 

fc=0 

where T(°)(n, FSi, FS2) = 1 and the higher order coefficients T('')(n, FSi, FS2) are given 
as polynomial expressions in A("')(n,FSi) and A(''')(n, FS2). 

For an arbitrary factorization scheme FS and arbitrary renormalization schemes RSi 
and RS2, there exists such a factorization scheme J^5(RSi, RS2, FS) that 

A(a;,M,FS,RSi) = A(x, M, J^cS(RSi, RS2, FS), RS2), (3.6) 

independently of x and the factorization scale M. The pair of {FS,RSi} thus defines 
the same singular factors that are absorbed into parton distribution functions as the pair 
of {J^5(RSi, RS2, FS), RS2}, and therefore both pairs define the same parton distribution 
functions and hard scattering cross-sections. This fact can be exploited for changing the 
renormalization scheme used for the factorization procedure because it allows to convert 
the simultaneous change of the factorization scheme and the renormalization scheme from 
{FSo, RSo} to {FS, RS} to changing the factorization scheme from FSq to J^5(RS, RSq, FS) 
for the fixed renormalization scheme RSq- 

The dependence of parton distribution functions on the factorization scale M is de- 
scribed by the evolution equations 

"^"^""d^'M '^^^ " ^' RS)D(n, M, FS, RS). (3.7) 

These equations are expressed in terms of Mellin moments. Converting them into x- 
space, we obtain the evolution equations in the form of (2.3). The splitting functions 
P(n, M, FS, RS) can be expanded in powers of a(M, RS) 

00 

P(n, M, FS, RS) = a^{M, RS) V^^\n, FS, RS). (3.8) 
fc=0 

Whereas the LO splitting functions p('')(n) are independent of the factorization scheme 
and the renormalization scheme, the higher order splitting functions 'P^^\n, FS, RS), k>l 
can be chosen at will and can be used for the specification of the appropriate factorization 
scheme.^ 



®If we specify factorization schemes via the corresponding higher order sphtting functions, then the 
complete specification of the factorization scheme requires also the specification of the corresponding renor- 
malization scheme because the relation between the splitting functions and the functions ^4'*'''' (a;), which 
define the factorization scheme, depends on the renormalization scheme. 



-9- 



The relation between T^^\n, FSi, FS2) for A; > 1 and the corresponding sphtting func- 
tions is given by the formula 

T('=)(n,FSi,FS2),P(°H^)l - A;6T('=)(n,FSi,FS2) = P^'^H", FSi, RS) - 

k-l 

-p('=)(n,FS2,RS) + ^ |p('=-')(n,FSi,RS)T(')(n,FSi,FS2) - 

1=1 

- T«(n, FSi, FS2)P('=~')(n, FS2, RS) + lbck-iiRS)T^^\n, FSi, FS2)} (3.9) 

where the definition of the coefficients q(RS), which are introduced in formula (2.1), is 
extended by co(RS) = 1 and ci(RS) = c. The preceding formula (3.9) holds for an arbi- 
trary renormalization scheme RS. If we specify factorization schemes via the corresponding 
splitting functions, then formula (3.9) allows us to determine the appropriate functions 
T(*'')(n, FSi, FS2), which are necessary for changing the factorization scheme (formula (3.9) 
forms a set of equations for T('^)(n, FSi, FS2), which can be solved iteratively).^*^ 

If higher order splitting functions are used for the specification of factorization schemes, 
then the formula for determining the factorization scheme J-"5(RSi, RS2, FS) has the form 

k 

P(^)(x, -F5(RSi, RS2, FS), RS2) = 4-/^(I^Si, RS2) P(')(x, FS, RSi), (3.10) 

1=0 

where the coefficients /ip^(RSi, RS2) are introduced in Appendix B.2. 
3.3 Hard scattering cross-sections 

This subsection contains an overview of formulae describing the dependence of hard scatter- 
ing cross-sections on the unphysical quantities associated with the factorization procedure. 
An arbitrary structure function F(x,Q^) is given as^^ 

F{n,Q^) = C(n,(5^M,FS,RS) D(n,M,FS,RS). (3.11) 

The coefficient functions C[n,Q^,M,FS,RS) can be expanded in powers of a{fi, RSe) 

00 

C(n,Q2,M,FS,RS) = J]]a'=(/i,RSE) C(*^)(n,Q^^,RSE,M,FS,RS). (3.12) 

fc=0 



^"To specify factorization schemes, we can also exploit functions T^*'^(n, FSi, FS2) — an arbitrary fac- 
torization scheme FS can be specified by the functions T^''\n, FS, FSo) where the factorization scheme FSo 
is some fixed and familiar factorization scheme. Formula (3.9) can then be used in the opposite way to 
determine the splitting functions in the factorization scheme FS (provided we know the splitting functions 
in the factorization scheme FSo). 

^^The coefficient functions form a row vector whereas the parton distribution functions are represented 
by a column vector, and therefore the multiplication on the right hand side of this equation yields a number 
(a matrix 1 x 1). 
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The dependence of the coefficient functions C'^^^(n, Q^, RSe, M, FS, RS) on the factor- 
ization scale M and the factorization scheme FS is governed by 

dCW(n,Q2 ,w,RS,M,FS,RS) ^ , 



dlnM 

1=0 



k-l-l 

(m+l) 



X E 5rr-i-i(^>/^'RS)p(-)(n,FS,RS), (3.13) 

m=0 

C('=)(n, /i, RS, M, FS, RS) = ^ C^^\n, Q^, fi, RS, M, FSq, RS) x 

1=0 
k-l 

X 9S_^(M,;U,RS)T(™)(n,FSo,FS), (3.14) 



where the renormalization scheme RSe, in which the renormalized couphng parameter 
used for expanding the coefficient functions is defined, is identical to the renormaHzation 
scheme RS, which is used for the factorization procedure. The coefficients (7;^'^^(/ii, /i2, RS) 
are introduced in Appendix B.2. The equivalence of the pair of {FS,RSi} to the pair of 
{-F5(RSi, RS2, FS), RS2} allows to convert changing the renormalization scheme used for 
the factorization procedure to changing the factorization scheme (for the fixed renormal- 
ization scheme), which is described by formula (3.14). Formulae (3.13) and (3.14) together 
with formula (B.7), which allows to change the renormalized coupling parameter used for 
expanding the coefficient functions, are thus sufficient for changing all unphysical param- 
eters associated with the renormalization and factorization procedure (even in the case if 
the renormalization scheme of the coupling parameter that is employed for expanding the 
coefficient functions is different from the renormalization scheme used for the factorization 
procedure) . 

Any inclusive cross-section cr{P) (in general differential) depending on observables P 
and describing a lepton-hadron collision is given as 

a{P) = V /" dx ai{x, P, M, FS, RS) A(x, M, FS, RS). (3.15) 
i Jo 

The hard scattering cross-section ai{x, P, M, FS, RS) can be expanded in powers of a{fi, RSe) 

00 

ai(x,P,M,FS,RS) = ^ a''+'=°(/i, RSe) fTf^(x, P, RSe, M, FS, RS) (3.16) 

fc=0 

where ko is a nonnegative integer. The formulae describing the dependence of the hard 

scattering cross-sections al''\x, P, RSe, M, FS, RS) on the factorization scale M and the 
factorization scheme FS read 

daf^(x,P,/x,RS,M, FS,RS) — /"^ — 
dlnM 

j ^1=0 



E / dy<^^af (xy,P,^,RSM,FS,RS) X 
' •'0 1 1=0 

-i-i ^ 
E 9ttym-iiM,n,RS)pjr\y,FS,RS)\ , (3.17) 

m=0 J 



k-l-l 

X 

m 
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1 ( ^ 

af )(x, P, fi, RS, M, FS, RS) = ^ / dy J ^ (xy, P, RS, M, FSq, RS) 

7 I «=o 



k-l ^ 

Y.gtljM,f,,RS)Tj:^\y,FSo,FS)\ . (3.18) 



X 

m=0 



Note that in the preceding formulae, the renormahzation scheme RSe is identical to the 
renormalization scheme RS. Formulae (3.17) and (3.18) represent an analogy of formulae 
(3.13) and (3.14) and together with formula (B.7) are sufficient for changing all unphysical 
quantities associated with the renormalization and factorization procedure. 

In the case of a hadron-hadron collision, any inclusive cross-section cr{P) depending 
on observables P can be expressed as 

-1 rl 

lo Jo 



(t{P) = Y,I I dxidx2aij(xi,X2,P,Mi,FSi,RSi,M2,FS2,RS2) X 

X D,/H^{xi,Mi,FSi,RSi)Dj/H,{x2,M2,FS2,RS2), (3.19) 



where the hard scattering cross-section (Tjj(xi, X2, -P, Mi, FSi, RSi, M2, FS2, RS2) can be 
expanded in powers of a(/u, RSe) 

c7ij(xi,X2,P,Mi,FSi,RSi,M2,FS2,RS2) = 

00 

= J2 I^Se) aff (xi , X2 , P, ^, RSe, Ml , FSi , RSi , M2 , FS2 , RS2) . (3.20) 

fe=0 

The dependence on the factorization scale Mi is determined by the formula 
dal.\xi,X2, P, RSi, Ml, FSi, RSi, M2, FS2, RS2) _ 



din Ml 



= -E / dyJ J^<TS^(xiy,x2,P,//,RSi,Mi,FSi,RSi,M2,FS2,RS2) X 

k-l-l •\ 

X J2 9tt-LliMl,^^,RSl)P;:^\y,FSl,RSl)\. (3.21) 

m=0 J 

The formula for changing the factorization scheme associated with hadron Hi (from FS^*^^ 
to FSi) has the form 

fjjf (xi , X2 , P, //, RSi , Ml , FSi , RSi , M2 , FS2 , RS2 ) = 



Yl [ dy|x;^i?(^iy>^2,P,/U,RSi,Mi,FsS°\RSi,M2,FS2,RS2) 



k-l \ 

X ^ gtljMuf^,RSi)T!^\y,FSf\FSi) . (3.22) 

m=0 J 



- 12 - 



Note that in the preceding formulae (3.21) and (3.22), the renormahzation scheme RSe is 
identified with the renormahzation scheme RSi. Analogous formulae hold for changing the 
unphysical quantities associated with hadron H2. The formula governing the dependence 
on the factorization scale M2 reads 

<lcj^^\xi,X2, P, RS2, Ml, FSi, RSi, M2, FS2, RS2) _ 



dlnM2 



Yl / dy<^ J^a5?(xi,X2y,P,/i,RS2,Mi,FSi,RSi,M2,FS2,RS2) X 

9tt-LiiM2,fi,RS2)P!:f{y,FS2,RS2)\. (3.23) 

m=0 J 



k-l-1 

X 

m= 



Changing the factorization scheme associated with hadron H2 (from FS2 to FS2) is de- 
termined by 

aff (xi , X2 , P, /X, RS2 , Ml , FSi , RSi , M2 , FS2 , RS2 ) = 



1 ( k 

= Y [ dyJ^cTj?(xi,X2y,P,//,RS2,Mi,FSi,RSi,M2,FS^°\RS2) X 

r -^0 I 1=0 

k-l ^ 

X Y ffS_™(M2,^,RS2)r(7)(y,FSf ,FS2) . (3.24) 

m=0 J 

Note that in the preceding formulae (3.23) and (3.24), the renormahzation scheme RSe 
is identified with the renormahzation scheme RS2. The above mentioned formulae (3.21), 
(3.22), (3.23), (3.24) together with formula (B.7) are sufficient for changing all unphysical 
quantities associated with the renormahzation and factorization procedure. 

3.4 Applicability of factorization schemes specified by splitting functions 

As it has already been mentioned, higher order splitting functions can be chosen at will 
and can be used for labeling factorization schemes. This subsection will be devoted to 
the question of practical applicability of factorization schemes that are specified by the 
corresponding splitting functions. 

To change the factorization scheme in which parton distribution functions and/or hard 
scattering cross-sections are defined, it is necessary to determine the appropriate functions 
T.^'^''(x, FSi, FS2). In the case when the factorization schemes are specified by the cor- 
responding higher order splitting functions, the necessary functions T;5,''^(x, FSi, FS2) are 
given as the solution of the system of equations represented by relation (3.9, D.30). The 
zeros of the denominators in formulae (C.3)-(C.9), which express the solution of the ap- 
propriate equations, can give rise to singularities in the Mellin moments T^j\n, FSi, FS2). 
The connection between the location of the singularitites of Mellin moments f{n) and the 
low X behaviour of the original function /(x), which is analysed in Appendix A, then imply 
that the zeros which are located sufficiently on the right (in the complex plane) can con- 
siderably influence the low x behaviour of the functions T-j\x, FSi, FS2). It can be proven 
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that for every real number ^, there exists such a real number kq that for every n > kq, the 
denominator given by (C.IO) has some zero point for Ren > ^, and therefore the influence 
of the zeros of the denominators on the low x behaviour of Tj^j^''(a;, FSi, FS2) cannot be 
ruled out, at least at higher orders (for large A;). 

Let us consider some factorization scheme FSq in which the low x behaviour of the ap- 
propriate parton distribution functions is fully determined by " physics" and is not affected 
by the choice of the finite parts A\j^\x), which define the factorization scheme. Singu- 
larities in the Mellin moments T^^^ {n,FS,FSo) induced by the zeros of the denominators 
may imply that the parton distribution functions in the factorization scheme FS have much 
larger values for low x than those in the factorization scheme FSq-^"^ If this occurs, then in 
the case of the factorization scheme FS, there must be an extensive mutual cancellation of 
those large values in the expressions for physical quantities which significantly depend on 
the low X region because the theoretical predictions for physical quantities have to be inde- 
pendent of the factorization scheme. The extensive mutual cancellation can cause problems 
in numerical calculations, and moreover, it is likely that the mutual cancellation is incom- 
plete at finite order calculations, which can result in unreliable theoretical predictions at 
the low X region. Hence, the range of the practical applicability of such a factorization 
scheme can be significantly restricted even though the corresponding splitting functions 
appear at first sight as reasonable. 

It is worth noting that if we specify factorization schemes using the appropriate finite 
parts A'Ij'^^x), then any unexpected restrictions of their practical applicability are ruled out 

because FS, MS) = Alf\x,FS), which follows from formulae (3.3), (D.5), (D.20) 

and (D.22). 

4 Situation at the next-to-leading order 

In the preceding section, we have presented formulae describing the dependence of parton 
distribution functions and hard scattering cross-sections on unphysical quantities associated 
with the factorization procedure. This section contains the compendium of these formulae 
and their consequences for NLO approximation, which consists in retaining only the first 
two terms in perturbative expansions. In this section, we limit ourselves only to the case 
when the renormalization scheme used for the factorization procedure is fixed and is the 
same as that in which the expansion parameter a(//, RS) is defined. The dependence on 
the renormalization scheme will thus not be written out explicitly. 

4.1 Changing the factorization scheme at NLO 

To change the factorization scheme at the NLO, it is necessary to determine the appro- 
priate functions Tj^ji'^^(x, FSi, FS2). According to formula (3.9, D.30), the Mellin moments 

^^The MS factorization scheme, in which the finite parts yl^*^"'(a;) are set equal to zero, should be an 
example of such a factorization scheme. 

^^Singularities in T^j \n, FS,FSo) that are induced by the zeros of the denominators may also cause 
undesirable behaviour of hard scattering cross-sections in the factorization scheme FS. 
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Tj^j\n, FSi, FS2) are given as the solution of the following equation 

Vi)(n,FSi,FS2),p(°Hn)l - 6TW(n, FSi, FS2) = pW(n, FSi) - pW(n, FS2). (4.1) 



The formula for expressing the functions T-j \x, FSi,FS2) in terms of the finite parts 
(x) then reads 

t(i)(x,FSi,FS2) = A(io)(x,FSi) - A(^°)(x,FS2), (4.2) 



which follows from relations (D.5) and (D.20). In NLO approximation, factorization 
schemes are thus fully specified by the corresponding finite parts ^^j'^^ (x) or NLO splitting 

functions P^^\x). If we use the functions A^^^\x) or T^^\x) for the specification of fac- 
torization schemes, then equation (4.1) allows us to determine the corresponding NLO 
splitting functions (provided we know them in some factorization scheme). 

4.2 Parton distribution functions 

The evolution equations in NLO approximation have the form 

"^^^dlnM^^ = (a(M)p(0)(n) + a2(M)p(i)(n,FS))D(n,M,FS). (4.3) 
The change of the factorization scheme is then described by the formula 

D(n,M,FS) = (^l + a(M)T(i)(n,FS,FSo))D(n,M,FSo), (4.4) 
which can be exponentiated to 

D(n,M,FS) =expfa(M)T(^)(n,FS,FSo))D(n,M,FSo). (4.5) 



An advantage of the exponentiated formula is the fact that the composition of FSq — )• 
FS and FS — )• FSq gives the original distribution functions, which does not hold for the 
unexponentiated formula. 

4.3 Hard scattering cross-sections 

In NLO approximation, the coefficient functions C(n, Q^, M, FS) are given as 

C(n,Q^M,FS) = C(°)(n,g2) +a(^)c(i)(n,Q^M,FS). (4.6) 

The LO coefficient functions C*-'^^(n, Q^) are independent of the unphysical quantities such 
as fjL, M and FS. The dependence of the NLO coefficient functions C^^^(n, Q^, M, FS) on 
the unphysical quantities is determined by 

C«(n,Q2^M,FS) = cW(n,Q2,Mo,FSo) + 

+ cW(n,Q2) (pW(n)ln^ + T«(n,FSo,FS)), (4.7) 



given factorization scheme FS can be specified by the functions T^j\x, FS, FSo) with the factorization 
scheme FSo to be some fixed and famihar factorization scheme (e.g. the MS factorization scheme). 
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which follows from relations (3.13, D.41), (3.14, D.42) and (B.7). 

Similarly, the hard scattering cross-section ai{x, P, M,FS) is expressed as 

ai(x,P,M,FS) = a'=«(;u)((jf^(x,P) + a(^)af^(x,P,/x,M,FS)) (4.8) 

where kg is some nonnegative integer. The NLO hard scattering cross-section a^^\x, P, n, M, FS) 
satisfies 

<Tf ^ (x, P, /i, M, FS) = a^P {x, P, fio, Mo, FSo) + hbln ^ (x, P) + 



^j\yaf\xy,P) (^P^W(y) In ^ + i;«(y, FSo, FS)) . (4.9) 



In the case of hadron-hadron collisions, we then have 



aij(xi,X2,P,Mi,FSi,Af2,FS2) 



= a'^'H/^) (<Tlf{xi,X2,P) + aifi) aff{xi,X2, P, /u. Mi, FSi, M2, FS2)) (4.10) 



with 



4'^(xi,X2,P,;U,Mi,FSi,M2,FS2) = 

= aff (xi , X2 , P, /xo , Mf ) , FSf) , M^ ^ , FS^^^ ) + A:o6 In ^ a^f (xi , X2 , P) + 

/io 

+ E/dy|^n (^iy'^2,P) (^P,f (y)ln^ + T(,^)(y,FSf),FSi)^ + 



+ ai°)(xi,X2y,P) (P,5)(y)ln^+r«(y,FSf ,FS2) ) }. (4.11) 



4.4 Applicability of factorization schemes specified by NLO splitting functions 

In Subsection 3.4, it has been shown that some splitting functions that appear at first 
sight as reasonable can correspond to factorization schemes which have some restrictions 
on their practical applicability. This subsection will therefore be devoted to a more detailed 
analysis of this fact at the NLO. 

A given factorization scheme FS can have some restrictions on its practical applicability 
if the Mellin moments Tj^j\n, FS,M.S) have some singularities that are located too much 
on the right in the complex plane. If we specify factorization schemes via NLO splitting 
functions, then such singularities can unexpectedly arise from the zeros of the denominators 
in formulae (C.3)-(C.9), which express the solution of equation (4.1). In the following, we 
limit ourselves only to the case when the matrix of the NLO splitting functions p(^)(x) has 
the same structure as that which corresponds to the MS factorization scheme, that is 
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In this case, the unexpected singularities in the corresponding Mellin moments Tj^j^ [n, FS, MS) 
can be induced only by the zeros of 

b' - W - - 8nfP(°)(n)p(°)(n). (4.13) 

A detailed analysis of the preceding expression for the number of quark flavours nf G 
{3,4,5} then shows that it has two simple roots in the half-plane Ren > 1. The approxi- 
mate numerical values of these roots are: 

n e {1.7329,4.6306} for nf = 3, 
n e {1.7995, 3.8458} for nf = 4, 

n G {1.9001,3.1798} for nf = 5. (4.14) 
Analysing formulae (C.3)-(C.9), which express the solution of equation (4.1), we find that 

r.(.i 

INTT.O snlittincr fnnrtinns / - - 



the Mellin moments T^j\n, FS, MS) do not have any pole at a simple root of (4.13) denoted 
by no if and only if the NLO splitting functions P^p{n) corresponding to the factorization 



scheme FS satisfy 



4°iK) (Pi?(no) - 4gK) - b) (p£)(no) - P^'Jim, MS)) + 



+ P(g (no) (PW (no) - P(S(^0) + b) (P« (no) - P^^ (no, MS)^ 
- 2p(g(no)P<g^(no) (P,«^(no) + ^^^^''(no) + 2nfP(i)^(no) - P^'^(no) 



P„«^(no, MS) - Pi,-^^(no, MS) - 2nfP„«^(no, MS) + P^ll(no, MS) ) = 0. (4.15) 



The preceding condition can be expressed in terms of the singlet splitting functions. 
At the NLO, the evolution of the quark singlet distribution function S(x) defined as 

^ix)=^{qiix) + qi{^)) (4.16) 
1=1 

and the gluon distribution function G{x) is described by the following system of coupled 
equations: 

^^1^ = {a{M)P^%{n) + a\M)P^'l{n)) S(n, M) + 
+ [a{M)P^%{n)+a\M)Pl^'l{n)) G{n,M), 

'^^^ = {a{M)p!,%)+a\M)P^'lin)) ^n^M) + 

+ {a{M)P§l{n) + a\M)pgl{n)) G{n, M) (4.17) 
where the singlet splitting functions are given by 

P(°^(x) = 2nfP(g(x), P^^lix) = 2nfP«(:r), 
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The condition (4.15) can then be rewritten in the form 

PSli^o) {pl,%{no) - pS(^o) - b) (P« (no) - P!,%o,m) + 
+ p('g(no) (P<g(no) - P<!S(no) + 6) (^^(^o) - 4q("o,MS)) - (4.19) 



- '^PhGi^o)P^Qi^o) [P^'^ino) - P'al,{n,) - P<^;^(no, MS) + P^;^(no, MS) j = 0, 

which means that this condition does not put any restriction on the choice of the non-singlet 
NLO sphtting functions. Hence, there should be no unexpected constraints on practical 
applicability in the non-singlet sector. 

5 Results of numerical analysis at NLO 

At the end of Section 2, we have shown that the ZERO factorization scheme could be useful 
for phenomenology and NLO Monte Carlo event generators. The important issue of its 
practical applicability at the NLO is the subject of this section. 



5.1 Changing the factorization scheme from MS to ZERO 

To convert hard scattering cross-sections and parton distribution functions from the stan- 
dard MS factorization scheme to the ZERO one, we need to know the functions T^j\x, MS, Zl 
which will be denoted as T^j \x) in the following. According to relation (4.1), the Mellin 
moments of these functions are given as the solution of 



T(i)(n),p(°)(n)] - 6TW(n) = P^^\n) (5.1) 

where P^^\x) denotes the MS NLO splitting functions [16, 17]. Using formulae (C.3)-(C.9), 
we get 



T^lin) = Tilin) = Tt\n) - J (^../^ ^^(n) + i^^)^(. 
TSkr^) = --A%)-2mT.^\n) (5.2) 



r^^H = Tilin) = T['\n), T§l{n) = Tg{n) = 4'\n), 
r«(n) = -^p(^)(n)-2nfr«(n) 



^^This is not surprising because the non-singlet case can be analysed in the same way as the general 
case in Appendix D with the only difference that all matrices are of dimension 1x1. Hence, formula (3.9, 
D.30) does not contain the commutator, and therefore no unexpected singularities can emerge in the Mellin 
moments i;^/^ (n, FSi, FS2). 

^®The functions Tj^^'^ {x, ZERO, MS), which are required for the conversion of parton distribution functions, 
are equal to MS, ZERO), which follows from formula (4.2). 
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where the Melhn moments T^^\n), T2^\n) and T^^\n) are given by (the dependence on 
n is not written out expHcitly in the foUowing formulae) 



(1) 



.(1) 



.(1) 



p(0) { p(0) _ p(0) 
^qG [-^qq ^GG 



p{i)V, pWy^ 2nfP(i)^ - P^^^ ] + 



(1) 



P, 



(1) 



Gq 



pSI {pS' - pSI + ^) {pS:^+ Pr+ - P^l) + 



+ (hP, 



(0) 



GG 



5p(0) 
<?<? 



6^ + 4nfPi°)pS,VS+4nf(P, 



qG Gq I ^ Gq 



,(0) 



p(0) P(0) _ p(0) , ,\ p(l) , p{0) /^p(O) _ p(0) 
-^gG \^qq ^GG ^ ] ^Gq ^ ^Gq \^qq ^GG 



h\ P. 



__ 2p(o)p(o) / p(i)y , p(i)^ , 2r7f P(i)^ - p(^) 

"^^qG^Gq yqq ^ ^qq ^ '^'n^qq ^GG 



gG 



(1) 



qG 



(5.3) 



(5.4) 



(5.5) 



and the denominator u is expressed as 



h b 



p(0) 
gg 



(0) 



GG 



Q„ p(0)p(0) 
^"f ^gG ^Gg 



(5.6) 



To obtain the functions T^p{x) in the x-space, it is necessary to determine the Melhn 

which has to be performed numerically. This was 



inversion of T^^\n), T2^'{n) and T, 
carried out for three and four (massless) quark flavours with the result 

T^'^\x) ^ dx-^ for x<0.1 



(1)^ 



n 



(5.7) 



where ^ = 4.63 for nf = 3 and ^ = 3.85 for Uf = 4. The values of the coefficients Cj are 
then such that the functions Tj^^\x) dominate over the NLO splitting functions P^j\x) for 
X < 0.1, which means that 



T^^l\x) ^ djx-'^ for x<0.1 



(5.8) 



This low x behaviour of the functions 



with the same value of ^ as in the case of T. 
Tj^j\x) is in agreement with the fact that the ZERO factorization scheme does not satisfy 

the condition (4.15). Since the rapid growth of the absolute value of T^^\x) occurs at 
X ~ 0.1, it is likely that the range of the practical applicability of the ZERO factorization 
scheme is restricted. 



5.2 The range of practical applicability of the ZERO factorization scheme 

The low X behaviour of the functions Tj^j\x, MS, ZERO) indicates that the ZERO factor- 
ization scheme has some restrictions on its practical applicability. In this subsection, we 
will investigate the practical applicability of the ZERO factorization scheme in the case of 
the structure function F2[x,Q'^Y 
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Figure 1. CoefRcient functions corresponding to the structure function F2[x,Q^) in the MS and 
ZERO factorization scheme for ni = 3. Note that because of the logarithmic scale, the coefficient 
functions in the right graph are plotted with the negative sign. 



A detailed numerical analysis of the ZERO factorization scheme was performed for 
three massless quark flavours. Using formula (4.5), the ZERO parton distribution func- 
tions at the factorization scale Mt = 100 GeV were calculated from the MS ones that 
corresponded to the MRST 1998 set [18]. The ZERO parton distribution functions at a 
general factorization scale M were then determined by the evolution from Mt to M in 
the ZERO factorization scheme. The theoretical predictions for i<2(x,(5^) were obtained 
by the numerically calculated Mellin inversion of the Mellin moments F2[n,Q'^^, which 
did not include any contributions concerning heavy quark flavours (only three (effectively) 
massless quark flavours were taken into account). Hence, the obtained theoretical predic- 
tions for -^2(3^) Q^) cannot be compared with experimental data, but are applicable for the 
assesment of the practical applicability of the ZERO factorization scheme, which is our 
aim. 

The right graph in Figure 1 shows that the quark singlet and gluon NLO coefficient 
function behave for low x in the same way as the functions T^^^^(x, MS, ZERO), and the 
same is true for the quark singlet and gluon distribution function, as can be seen from 
Figure 2. The low x behaviour of these distribution functions implies that if some physical 
quantity depends on the values of these distribution functions in the region x < 0.1, 
then obtaining a reasonable theoretical prediction for this quantity requires a considerable 
mutual cancellation of large values in the appropriate formula. It is likely that the sufficient 
cancellation occurs, if ever, only for some choices of the renormalization and factorization 
scale, which means that most of the choices result in unreliable theoretical predictions. 
Moreover, the sufficient cancellation causes complications in numerical computations, which 
can be hardly soluble. 

Both undesirable facts manifest themselves in Figure 3. The upper graphs show a 
theoretical prediction for F2[x,Q'^), which is unreliable for x < 0.1. In this case, the the- 
oretically predicted F2[x,Q^) rapidly grows with decreasing x in the region x < 0.1. The 
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Figure 2. Quark singlet and gluon distribution function in the MS and ZERO factorization scheme 
for n{ — 3. Note that the gluon distribution functions in the lower right graph are plotted in their 
absolute value — whereas the MS distribution is positive, the ZERO distribution is negative (see 
the lower left graph). 

lower left graph indicates that there are some theoretical predictions for i<2(a;,(5^) which 
rapidly fall with decreasing x. Hence, the range of the theoretical predictions for F2[x, Q^) 
is very large for x < 0.1, and therefore the ZERO factorization scheme for nf = 3 has 
practically no predictive power in the region x < 0.1. The range of the theoretical pre- 
dictions for F2[x,Q'^^ should include reasonable values, which means that for a given xq, 
there exist such choices of the renormalization and factorization scale that the correspond- 
ing theoretical predictions for the value of F2[xo,Q'^^ are reasonable (but this does not 
mean that there must be such a choice of the renormalization and factorization scale that 
the corresponding theoretical prediction for F2[x,Q'^^ is reasonable for all x). However, 
obtaining the reasonable theoretical predictions in numerical calculations is difficult, which 
can be seen in the lower right graph, where the used precision of the numerical compu- 
tation is insufficient for x < 0.02. The comparison of the lower graphs shows that the 
reasonable theoretical predictions in the region a; < 0.1 are very sensitive to the choice of 
the renormalization and factorization scale, which is a consequence of the sensitivity of the 
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Figure 3. Theoretical predictions for F2(^x,Q^) in the MS and ZERO factorization scheme for 
nf — 3. 




Figure 4. Quark non-singlet distribution function and the non-singlet contribution -F'2,ns (a^, Q^) 
to -^2(2;, Q^) in the MS and ZERO factorization scheme for rif = 3. 
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Figure 5. Quark singlet and gluon distribution function in the MS and ZERO factorization scheme 
for fif — 5. Note that the gluon distribution functions in the lower right graph are plotted in their 
absolute value — whereas the MS distribution is positive, the ZERO distribution is negative (see 
the lower left graph). 



extent of the mutual cancellation of large values to this choice. 

According to Subsection 4.4, the ZERO factorization scheme should have no restric- 
tions on its practical applicability in the non-singlet case, which is demonstrated in Figure 4. 

The functions r/^^^(3;,MS, ZERO) were not calculated in the 2;-space for five (massless) 
quark flavours. Since the ZERO factorization scheme for rif = 5 does not satisfy the 
condition (4.15), the functions T-j \x, MS, ZERO) should behave for low x in the following 
way (see (4.14)): 

Tij\x,MS, ZERO) « djx-^, ^ = 3.18, (5.9) 

which indicates that the ZERO factorization scheme has some constraints on its practical 
applicability, similarly as in the case of rif = 3 and nf = 4. The ZERO parton distribution 
functions and the theoretical predictions for F2[x, Q^) were calculated in a similar way^'' as 

^^Seeing that the Mellin moments of the parton distribution functions and the structure function 
F2{x,Q^) tend to zero sufficiently fast for large n, the numerical computation of their Mellin inversion 
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Figure 6. Theoretical predictions for F2(^x,Q^) in the MS and ZERO factorization scheme for 
nf — 5. 

in the case of nf = 3 (the only difference was the fact that the value of the factorization scale 
Mt was lO'^ GeV and the MS parton distribution functions corresponded to the MSTW 
2008 set [19]). The ZERO parton distribution functions displayed in Figure 5 behave for 
low j; in a similar way as in the case of rif = 3, which is displayed in Figure 2. The main 
difference is the rate of the growth, which is in accordance with relation (5.9) and is thus 
slower than in the case of nf = 3. The comparison of Figures 2 and 5 also shows that 
in the case of nf = 5, the rapid growth occurs at a bit lower value of x. However, the 
difference is not significant. Hence, it is not surprising that the ZERO factorization scheme 
for nf = 5 gives unreliable theoretical predictions for i<2(x,Q^) in the region of low x, as 
can be seen in Figure 6. The comparison of the upper and lower graphs in Figure 6 shows 
that the region of the practical applicability increases slowly with increasing energy, which 
also holds for nf = 3. 

For the number of massless quark flavours from three to five, which represents the 
relevant values for QCD phenomenology, the ZERO factorization scheme is fully applica- 

is quite easy, which does not hold for the Mellin inversion of 7;^''(n,MS,ZER0). 
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ble for the physical quantities that do not depend on the values of the quark singlet and 
gluon distribution function in the region x < xq where xq = 0.1. However, in the case 
of the physical quantities which depend on the values of these distribution functions in 
the region x < xq, the ZERO factorization scheme gives unreliable theoretical predictions 
and has little predictive power. Even for these physical quantities, obtaining reasonable 
theoretical predictions is possible for some special choices of the renormalization and fac- 
torization scale, but requires a considerable mutual cancellation of large values, which 
causes difficulties in numerical calculations. In some cases, the value of the bound xq can 
be lowered, however, it is practically impossible to lower the value of xq in such a way 
that the ZERO factorization scheme would be applicable in the full range that is used in 
QCD phenomenology because the extent of the undesirable features rapidly grows with 
decreasing xq. 

As it has already been mentioned, the ZERO factorization scheme puts all NLO cor- 
rections into hard scattering cross-sections and is therefore a certain opposite to the DIS 
factorization scheme, in which all NLO corrections (to i^2(x,(5^)) are included in the cor- 
responding NLO splitting functions. The restricted practical applicability of the ZERO 
factorization scheme is thus surprising because the DIS factorization scheme can be ap- 
plied without any restrictions. 

There are other reasons why the restricted practical applicability of the ZERO factor- 
ization scheme is unexpected and not obvious. As it has already been shown, the ZERO 
factorization scheme has no restrictions of its practical applicability in the non-singlet case. 
In addition, the condition of practical applicability (eq. (4.15)) is nontrivial and its conse- 
quences depend on the number of massless quark flavours (e.g. the growth of the ZERO 
quark singlet distribution function in the low x region is a result of the fact that the ZERO 
factorization scheme does not satisfy the condition of practical applicability and the rate 
of this growth depends nontrivially on the number of massless quark flavours (see Figures 
2 and 5)). 

6 Summary and conclusion 

In this text, we have analysed the freedom associated with factorization in massless per- 
turbative QCD, with the emphasis on its quantification. We have derived the formulae 
that are useful for changing the unphysical quantities associated with the renormalization 
and factorization procedure and we have shown that factorization schemes can be specified 
using the higher order splitting functions, which can be chosen at will. This allows us to 
introduce the so called ZERO factorization scheme, which is defined by setting the higher 
order splitting functions equal to zero, which means that the ZERO factorization scheme 
represents a certain opposite to the more familiar DIS factorization scheme. The poten- 
tial exploitation of the ZERO factorization scheme for the construction of consistent NLO 
Monte Carlo event generators in which initial state parton showers can be taken formally 
at the LO has provided the main motivation for this study. The related questions have 
been discussed in more detail at the end of Section 2. 
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The discussion of the practical apphcabihty of general factorization schemes has shown 
that the factorization schemes specified by splitting functions can have unexpected restric- 
tions on its practical applicability. Even if splitting functions appear at first sight as 
reasonable, they can specify such a factorization scheme that the values of the appropriate 
parton distribution functions are significantly larger in the low x region than in the case 
of the standard MS factorization scheme. The cancellation of these large values, which has 
to occur in the expressions for physical quantities, can then result in unreliable theoretical 
predictions, little predictive power and difficulties in numerical calculations, which restricts 
the practical applicability of such a factorization scheme. 

The question of the practical applicability of general factorization schemes has been 
studied in more detail at the NLO. We have found the condition (eq. (4.15)) that must 
be satisfied by NLO splitting functions in order that any unexpected restrictions on the 
practical applicability of the corresponding factorization scheme are ruled out. Unfortu- 
nately, this condition is not satisfied for the ZERO factorization scheme. The numerical 
analysis of the ZERO factorization scheme at the NLO has then shown that its application 
is reliable for the physical quantities which do not depend on the values of the quark singlet 
and gluon distribution function in the region x < 0.1. However, if we apply the ZERO 
factorization scheme for the other physical quantities (i.e. those which depend on the values 
of these distribution functions in the region x < 0.1), then the undesirable features men- 
tioned in the preceding paragraph can occur. The restricted practical applicability of the 
ZERO factorization scheme has been illustrated on the structure function F2[x,Q'^^ in the 
low X region. However, even for the physical quantities for which the ZERO factorization 
scheme is unreliable, there should exist some special choices of the renormalization and 
factorization scale that result in reasonable theoretical predictions, but the exploitation 
of this fact in Monte Carlo event generators is practically impossible because obtaining 
reasonable theoretical predictions in this case requires a considerable mutual cancellation 
of large values. 

Although the ZERO factorization scheme has some restrictions of its practical appli- 
cability, its phenomenological exploitation still makes sense because it can be applied for 
the description of the production of heavy objects, which is important in searches for new 
physics. 

A Mellin transform and related formulae 

The Mellin transform of a function f{x) defined on the interval (0, 1) is given as 



where n is in general a complex number. The preceding integral, if exists at least for some 
real value of n, defines a function of complex variable n that is holomorphic at least in 
some right half plane. Using analytical continuation, the Mellin moments /(n) (i.e. the 
Mellin transform of the function /(x))^^ can be defined even for some values of n for which 

^*In this text, the same symbol is used for both the function and its MelHn moments. The distinction 
between them is made by the argument — using the argument n always refers to the Mellin moments. 
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the integral does not exist (however, this does not mean that the Mehin moments can be 
defined for all n). The Mellin moments /(n) uniquely determine the original function /(x), 
and therefore the Mellin transform can be inverted. The Mellin transform and its inversion 
are linear. The formula for the inverse Mellin transform reads 



1 /"5+ioo 

fix) = — / /(n)x-"dn (A.2) 

' ^— ioo 
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where ^ is a real parameter chosen in such a way that the integration contour is located to 
the right of all singularities of /(n) in the plane of complex n. 

The location of the singularities of the Mellin moments fin) is related to the behaviour 
of the function fix) in the vicinity of x = 0. If the Mellin moments fin) have some 
singularity for Re n > and if there exists a nonnegative integer k such that the function 
gix) = (1 — x)^ fix) is bounded on some interval [e, 1], then there is no real number C such 
that ^ 

\fix)\<c{^ forxG(0,e) (A.3) 

because if the preceding relation held for some real number C, then the function gix) would 
have to satisfy 



\gix)\<K[-^J (A.4) 

for some real number K and for all x from the interval (0, 1]. This constraint would then 
imply the following bound 

\gin)\< C \x''-^\\gix)\Ax < K C x^'^^'-^-^dx, (A.5) 
Jo Jo 

which would mean that the Mellin moments gin) could not have any singularity for Ren > 
but this would be in contradiction with the assumption that the Mellin moments fin) 
have some singularity for Re n > ^ because 

gin) = f ^(-1)^- x^fix) dx = ^(-1^ /(n + j). (A.6) 

Some relations concerning factorization (e.g. the relation between dressed and bare 
parton distribution functions) are expressed as a convolution given by 

(/ * 9)ix) = C f'dydz fiy)giz)6iyz - x) = (A.7) 



JO 

''^fi^)s{y)= t'^msC^)- (A.8) 

y \yJ Jx y \yJ 

From formula (A.7), we easily get 

if*9)in) = fin)gin), (A.9) 
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which means that the Mehin transform converts the convolution (A. 7) to ordinary multi- 
phcation of Melhn moments. 

Using formula (A. 7), we find 

/ d.xf{x){g*h){x)= / / dydz f{yz)g{y)h{z) = 
Jo Jo Jo 

1 .1 

dx {f (g) g){x)h{x) = / dx {f (g) h){x)g{x) (A.IO) 
Jo 

where the symbol (8> is defined as 

if'gg){x) = [ dyf{xy)g{y) (A.ll) 
Jo 

and satisfies 

[{f ^ g) ^ h]{x) = I dz{f(gg){xz)h{z)= [ [ dydz f{xyz)g{y)h{z) = 
Jo Jo Jo 

= [ [ dydzg{y)h{z) [ dU{xOHC-yz) = 
Jo Jo Jo 

dUixO [ I dydzg{y)h{z)5{i-yz) = 
Jo Jo 

di f{xO {g * hm = [/ (5 * h)]{x). (A.12) 



Formulae (A.IO) and (A.12) are useful for manipulations with expressions for cross-sections. 

Using relations (A.l), (A. 7) and (A.ll), we easily get the following properties of the 
function f{x) = 5{1 — x): 

/(n) = l, f*g = g*f = g, g®f = g (A.13) 

where g is an arbitrary function. 



B Some important formulae for the QCD coupling parameter a(/i, RS) 
B.l Changing expansion parameter 

Within the framework of perturbative QCD, the QCD coupling parameter a(/i, RS) is a 
function of the renormalization scale ^ and the renormalization scheme RS. Hence, it is 
useful to be familiar with relations between coupling parameters corresponding to different 
arguments. 

Consider an expansion parameter a{P) depending on a set of parameters denoted by P. 
Let us assume that the relation between expansion parameters corresponding to different 
sets of parameters Pi and P2 has the form of power series: 

00 

a{Pi) = Y,cf\Pi,P2)o'+\P2). (B.l) 
1=0 
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This relation can then be generahzed as follows: 

oo 

aHPi) = Y.4''\Pi,P2)a'+HP2) (B.2) 



(k) 

c 

1=0 

where fc is a nonnegative integer and 



cf\Pi,P2) = 6io, (B.3) 



-I 

I 



4'^'\p^,P2) = Y,4-JPi^P2)cgHPi,P2) for k>l. (B.4) 

m=0 

Formula (B.2) is useful for changing an expansion parameter of power expansions. Consider 
some quantity F which can be expressed as an expansion in powers of a(Po): 

oo 

^ = ^^(fc)(Pp)„fc+fco(pp), (B.5) 

k=0 

Using formula (B.2), we can then obtain an expansion of F in powers of an arbitrary 
expansion parameter a{P)- 

oo 

F = Y, F^''\P) a^+^''{P) (B.6) 

k=0 

where 



k 

= ^F«(Po)4'+f°)(Po,P). (B.7) 

1=0 

B.2 The QCD coupling parameter in d = 4 — 2e space-time dimensions 

Within the framework of dimensional regularization in d = 4 — 2e space-time dimensions, 
equation (2.1), which describes the dependence of the QCD coupling parameter a(/x, RS) 
on the renormalization scale /x, has the following form 



^^^^ = f;A(RS)a'+i(^,RS) (B.^ 



din// 

where /3o(RS) = — 2e and all higher order coefficients /3/(RS), / > 1 are linear functions of 
£ and depend on the renormalization scheme RS (in the case of the coefficients /3i(RS) and 
/?2(RS), only the part proportional to e is dependent on the renormalization scheme RS). 

The renormalization scale of the QCD coupling parameter can be changed using the 
formula 

oo 

which can be derived from (B.8). The coefficients^^ g^j^\^, depend on e and are 

finite in the limit e — )• 0. The change of the renormalization scheme of the QCD coupling 



l{^i,m) = Y,gf\^i,^l^,RS)a'+\^JiQ,RS), (b.9) 



^^These coefficients can also be written as gi (^,RS) because they depend only on the ratio of ^ and 
Mo- 
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parameter is then described by the formula^*^ 

oo 

a(^,RSi) = ^/i[^^(RSi,RS2)a'+H/i,RS2), /i[,^^(RSi, RS2) = 1. (B.IO) 

1=0 

Contrary to the coefficients g fj, fiQ, RS), the coefficients /ip^(RSi, RS2) are independent 
of e. Both formulae (B.9) and (B.IO) have the form of (B.l) and can thus be generahzed 



to 

00 



a^X/i,RS) = j;gf)(^,Mo,RS)a'+'=(/io,RS), (B.ll) 



1=0 

00 

'=(/i,RSi) = J];/ip^(RSi,RS2)a'+^XM,RS2), (B.12) 



a _ 

1=0 

where the coefficients gl''\^, fio,KS) and /ip^(RSi, RS2) are determined by relations (B.3) 
and (B.4). 

C Solution of equations (D.26) and (3.9, D.30) 

This appendix contains formulae for the solution of the matrix equation 

X(n),p(°)(n)j -KX(n) = Y(n), (C.l) 

with respect to X(n). The matrix p('')(n) represents the LO splitting functions and the 
matrix Y(n) is given. All three matrices in the preceding equation are square matrices 
with indices corresponding to parton species: qi, qi and G, and therefore their dimension 
is 2nf + 1 where nf denotes the number of quark flavours. For the following, it is useful to 
introduce a unified labeling for quarks and antiquarks: Qi = qi, Qi+Uf = qi- The structure 
of the matrix P^'^\n) can then be expressed as follows: 

For better readability of the formulae expressing the solution of equation (C.l), it is 
convenient not to write out explicitly the dependence on n. The formulae then read 

(0) p(0) / ,^2nt ^ V \ ,^ p(0) 



X,.,^ - X - -Y,^,^ + ^ ^^^^^ _ ^^^^ ^ ^^^^^^^^^^ + 



p^dp'cl Efil yq.q, - y,,) - <j [ycQ, - yc, 



^ { yQ.G - y,G) + p$ yq.q. - y,, 

'^'^GG '^^ll ^ ^ '^^^^qG ^Gq 



(C.3) 



Xq.g = X,g + ^(0) ^p(o) ,.2 , o.„p(0)p(0) ' (^-^^ 

■^"Note that the coefBcients /i^' (RSi, RS2) do not depend on ^. 
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Xgg 



1 



k{Ygq,-Yg,)-P, 



(0) 
Gq 



E2nf 
k=i ^QkQi 



(0) 



Ygg - 2niX 



(0) 



GG 



K^ + 2n{P'"'P, 



(0) d(0) 



qG-" Gq 



where 



X 



1 r 



>(0) (p{0) 
Gq y^qq 



Gq 
3(0) 

Gq 



P, 



(0) 



GG 



2p(0)p{0) ^y^^ 



Ygg) 



X. 



qG 



^Gq 



1 r 



- - {^yqG + ^G iYqq " ^^g)) + 



1 r 

1/ 



p, 



(0) 



GG 



p{0) 



99 



^^Gq - P^gI iYqq " ^Gg) ] + 



+ 4n<)(p£)y,c + pi°VG. 



and the denominator v is given as 



K K 



p(0) 



p(o)y 

^GGJ 



Q„ p(0)p(0) 
^^f PgG ^Gq 



The symbols Fgq, and Ygq are defined as follows: 

^ 2nf ^ 2rif ^ 2nf 



k,l=l 



2ni 



k=l 



2ni 



GQk 



k=l 



(C.5) 
(C.6) 



(C.7) 



(C.8) 



(C.9) 



(C.IO) 



(C.ll) 



If all denominators in expressions (C.3)-(C.9) are nonzero, then the solution of equation 
(C.l) exists and is unique for an arbitrary right hand side Y. 

It can be shown that for every nonzero k, there exists some right half plane in which all 
denominators in expressions (C.3)-(C.9) are nonzero. Hence, if the right hand side Y(n) 
of (C.l) is holomorphic in some right half plane, then there exists some right half plane in 
which the solution X(n) of equation (C.l) is unique and holomorphic. 



D Analysis of the freedom associated with factorization in massless QCD 

This appendix contains a detailed analysis of the freedom associated with the factorization 
procedure in massless perturbative QCD. Important results of this appendix are summa- 
rized in Subsections 3.2 and 3.3. 

To investigate the ambiguity associated with the definition of dressed parton distribu- 
tion functions, it is useful to express the relevant formulae in terms of Mellin moments, 
which are defined in Appendix A. Converting relation (3.1) into the space of Mellin mo- 
ments, we get ordinary matrix multiplication of moments 

D(n, M, FS, RS) = A(n, M, FS, RS)D(n) (D.l) 
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and the expansion of A(n, M, FS,RS) is expressed as follows 

oo 

A(n,M,FS,RS) = ^a'=(M,RS)A('=)(n,FS,RS), A(°)(n, FS, RS) = 1. (D.2) 

fc=0 

The form of relation (D.l) allows us to obtain expression for the bare parton distribution 
functions D(n): 

D(n) = B(n, M, FS, RS)D(n, M, FS, RS) (D.3) 

where B(n, M, FS, RS) is the matrix inversion of A(n, M, FS, RS), which can be expanded 
in powers of a(M, RS): 

oo 

B(n, M, FS, RS) = a^M, RS)B('=)(n, FS, RS) (D.4) 

where the coefficients B'^^''(n, FS, RS) are determined by the following recurrence relation 

jfc-i 

B(*^)(n,FS,RS) = - J]] A('=-')(n,FS,RS)B(')(n,FS,RS), B(°)(n, FS, RS) = 1, (D.5) 

1=0 

which implies that 'B^^\n, FS, RS) is given as a sum of terms of the form 

cA('i)(n, FS, RS) A('2)(n, FS, RS) • • • A('-)(n, FS, RS) (D.6) 

m 

where the coefficient c is an integer and li = k. This means that B^'^) (n, FS, RS) depends 

i=l 

on A(')(n,FS,RS) for 1 < / < /t and is independent of A(')(n, FS, RS) for / > k. Relation 
(D.5) and its consequence are useful for analyzing formulae (D.14) and (D.20). 

The following three subsections will be devoted to changing the unphysical quantities 
on which the dressed parton distribution functions depend. The obtained results will then 
allow us to investigate the dependence of the splitting functions and the hard scattering 
cross-sections on the unphysical quantities associated with the factorization procedure. 

D.l Changing the renormalization scheme 

If we want to change the renormalization scheme used for the factorization procedure from 
RSi to RS2, then we should convert the expansion parameter a(M, RSi) to a(M, RS2). 
Hence, let us express Aij{x, M, FS, RSi) as an expansion in powers of a(M, RS2): 

A,j{x, M,FS,RSi) = 6ij6{l-x) + 

00 k 

+ a^M, RS2) 4-z(RSi, RS2) Af^{x, FS, RSi), (D.7) 

k=l 1=1 

where we have used formula (B.12). Seeing that the coefficients /i['^^(RSi, RS2) are inde- 
pendent of £, the preceding expansion has the form that is consistent with formulae (3.2) 
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and (3.3), and therefore there must exist such a factorization scheme J-'5(RSi, RS2, FS) 
that 

Aij{x, M,FS,RSi) = ^,j(x,M,J^5(RSi,RS2,FS),RS2). (D.8) 

The pair of {FS, RSi} thus defines, independently of the factorization scale M, the same sin- 
gular factors that are absorbed into parton distribution functions as the pair of {J^5(RSi, RS2, FS), RS2}, 
and therefore both pairs are equivalent in the sense of the definition of dressed parton dis- 
tribution functions. Changing the renormalization scheme from RSq to RS for the fixed 
factorization scheme FS and factorization scale M can thus be converted to changing the 
factorization scheme FS to the factorization scheme J-"5(RS, RSq, FS) for the fixed renor- 
malization scheme RSq and factorization scale M. The change of the factorization scheme 
will be discussed in Subsection D.3. 

In order to be able to change the renormalization scheme in the above described man- 
ner, we have to know some prescription for the determination of the factorization scheme 
-F5(RSi, RS2, FS). Such a prescription can be obtained as follows. Comparing the expan- 
sion (3.2) of the right hand side of (D.8) with the expansion (D.7) of the left hand side of 
(D.8), we find 

k 

4j^^(x, J-5(RSi,RS2,FS),RS2) = 4-i(i^Si, RS2) ^^^(x, FS, RSi), k>l. (D.9) 

1=1 

Using formula (3.3) and exploiting the fact that the coefficients /i['^''(RSi, RS2) are inde- 
pendent of e, we get 



4^^"'^(a;,^cS(RSi,RS2,FS),RS2) = ^ /ifl;(RSi, RS2) aS5'"^(x, FS, RSi) (D.IO) 



T3G_ T5G. T?Q\ ^^C.^— ^^^W T3Q.^ /I (''") / 

1=1 

where k > 1 and m > 0. The desired prescription is then represented by the preceding 
equation (D.IO) for m = 0. 

D.2 Evolution equations 

The evolution equations describe the dependence of dressed parton distribution functions 
on the factorization scale. From formula (D.l), we get 

dD(n,M,FS,RS) dA(n, M, FS, RS) t,q rq^ 

= B(n, M, FS, RS)D(n, M, FS, RS), (D.ll) 

where we have used formula (D.3) to express the bare parton distribution functions D(n) in 
terms of the dressed parton distribution functions D(n, M, FS, RS). Expanding in powers 
of a(M, RS), using equation (B.8) and after some manipulations, we obtain 

^^^^^^j^|^^^^^B(n, M, FS, RS) = a(M, RS)P(n, M, FS, RS) (D.12) 

where the splitting functions P(n, M, FS, RS) have an expansion in powers of a(M, RS) 

00 

P(n,M,FS,RS) = J]]a*^(M,RS)p(*^)(n,FS,RS) (D.13) 

k=0 
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with the coefficients given by the formula 

p('^)(n,FS,RS) = ^ I ^(p+l)A-p(RS)A(P+i)(n,FS,RS) j B('=-')(n, FS, RS). (D.14) 
1=0 \p=0 J 

The evolution equations expressed in terms of Mellin moments thus read 

dD(n M,FS,RS) ^ RS)p(n, M, FS, RS)D(n, M, FS, RS). (D.15) 
dmM 

Converting the preceding equation into x-space, we obtain the evolution equations in the 
form of (2.3). 

The splitting functions P^'^^(n, FS, RS), determined by (D.14), must be finite in the 
limit e — 7- 0. A detailed analysis of this condition shows that 

A(fcO(n,FS,RS) = for / > A; > 1, 

A(^')(n, FS, RS) for A; > / > 2 is uniquely determined by (D.16) 
A(™o)(n,FS) and A(™i)(n, FS, RS) for 1 < m < /t - 1. 

Since the functions A^j'^\x) can be choosen at will, all nontrivial information concerning 
the properties of the theory has to be contained in the functions A^-j^\x, FS, RS). 

D.3 Changing the factorization scheme 

The formula for changing the factorization scheme can be obtained from relations (D.l) 
and (D.3): 

D(n, M, FSi, RS) = T(n, M, FSi, FS2, RS)D(n, M, FS2, RS) (D.17) 

where 

T(n, M, FSi, FS2, RS) = A(n, M, FSi, RS)B(n, M, FS2, RS). (D.18) 

Expanding the preceding relation in powers of a(M, RS), we get 

00 

T(n, M, FSi, FS2, RS) = ^ a'=(M, RS) T^''\n, FSi, FS2, RS) (D.19) 
where the coefficients of the expansion are determined by the formula 

k 

T(*^)(n,FSi,FS2,RS) = ^ A('=-')(n, FSi, RS)B(')(n, FS2, RS). (D.20) 

/=o 

Using relation (D.5), we find that T(o)(n, FSi, FS2, RS) = 1 and T('=)(n, FSi, FS2, RS) 
for /c > 1 is given as a polynomial expression in AW(n,FSi,RS) and A«(n,FS2,RS) for 
1 < I < k with coefficients independent of e. Seeing that T('=)(n,FSi,FS2,RS) must be 
finite in the limit e — )• 0, we can do the following replacement^^ 

A(')(n,FSi,RS) ^ A('°)(n,FSi), A(')(n, FS2, RS) ^ A('°)(n, FS2) (D.21) 



■^^The functions Aij {x), which define the factorization scheme, do not depend on the renormahzation 
scheme. 
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in the expression for T('=)(n, FSi, FS2, RS) because the singular parts of A(')(n, FSi, RS) and 
AW(n,FS2,RS) cannot contribute to the finite part of T('=)(n, FSi, FS2, RS) (see formula 
(3.3)). Hence, T('^)(n, FSi, FS2, RS) does not depend on the renormalization scheme RS. 

The condition that T(^)(?7-, FSi, FS2) is finite in the limit e — t- allows us to determine 
the dependence of A('^)(n, FS, RS) on the factorization scheme. Consider the so called 
minimal subtraction (MS) factorization scheme, which is defined by setting the functions 
A'Ij^\x) equal to zero. From relations (3.3) and (D.5), we get 

B^'') (n, MS, RS) = - - A^'^^) (n, MS, RS) + terms containing at least \ . (D.22) 

The requirement of the absence of a term proportional to 1/e on the right hand side of 
(D.20) then implies 

k-l 

A(^^) (n, FS, RS) = A(") (n, MS, RS) + ^ A^'^-'-O) (n, FS) A^'^) (n, MS, RS) . (D.23) 

1=1 

This formula determines the dependence of the functions A^j^\x, FS, RS) on the arbitrary 

functions A^-j^^ (x), which specify the factorization scheme. The dependence of the functions 

A'^j^\x, MS, RS) on the renormalization scheme is then described by the following formula 
resulting from (D.IO): 

fc 

A\f\x, MS, RS) = J2 4-/(RSo, RS) Afl\x, MS, RSq), (D.24) 
1=1 

where we have used the fact that J^5(RSo, RS, MS) = MS, which follows from equation 
(D.IO) for m = 0. The obtained relations (D.16), (D.23) and (D.24) thus imply that if we 
know the functions A^j^^(x, FS, RS) for some factorization scheme and some renormaliza- 
tion scheme, then we are able to determine the complete singular factors Aij{x, M, FS, RS) 
for an arbitrary factorization scale, factorization scheme and renormalization scheme. 

Since the pair of {FS,RS} is equivalent to the pair of {J-'5(RS, RSq, FS), RSq}, we 
can convert the simultaneous change of the factorization scheme and the renormalization 
scheme from {FSo,RSo} to {FS,RS} to changing the factorization scheme from FSq to 
-F5(RS, RSo, FS) for the fixed renormalization scheme RSq. 

D.4 Specification of factorization scheme via splitting functions 

The relations presented in the preceding subsections allow us to investigate the dependence 
of the splitting functions on the factorization scheme, which is one of the main aims of this 
text. This dependence interests us only in the limit e — )• (in four space-time dimensions). 
From formula (D.14), we get that the LO splitting functions for e = are given as 

p(°Vx,FS,RS) = lim f-2eA(i)(x, FS, RS)) = -2A(ii)(x, FS, RS). (D.25) 

Relations (D.23), (D.24) and (B.IO) then imply that A(ii)(a;, FS, RS) is independent of 
the factorization scheme and the renormalization scheme, and therefore the LO splitting 
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functions P*-'^^(a;) are unique (for e = 0). However, this is not true for the higher order 
sphtting functions. Using (D.5), (D.14), (D.23), (D.25) and the properties of the coeffi- 
cients /3fc(RS) mentioned in Appendix B.2, we find that the higher order sphtting functions 
p(^)(n, FS, RS), A; > 1 for e = are expressed as 

p('=)(n,FS,RS) = [A('=°)(n,FS),p(°)(n)j - fc^A^^^^^K FS) + Rfc(n, FS, RS) (D.26) 

where Rfc(n, FS,RS) is a polynomial expression in A('i)(n,MS,RS) for 1 < Z < A; + 1 and 
A^^^\n, FS) for 1 < I < k — 1 with the coefficients depending on the renormalization scheme 
RS. The preceding relation (D.26) and the relations presented in Appendix C imply that 
for any set of arbitrarily chosen splitting functions p(')(n), 1 < I < k where k is arbitrary, 
there exists exactly one set of the corresponding A("^)(n), 1 < I < k. The higher order 
splitting functions P('^)(x, FS, RS), k > 1 for e = can thus be chosen at will and can be 
used for labeling factorization schemes. 

The specification of factorization schemes via splitting functions can be useful for phe- 
nomenology. In order to simplify the application of this kind of specification, it is necessary 
to derive formulae for expressing T('^)(x, FSi, FS2) and -F5(RSi, RS2, FS) in terms of split- 
ting functions. The relation between T('^)(2;, FSi, FS2) and the corresponding splitting 
functions can be obtained as follows. Let us start with the relation 

A(n, M, FSi, RS) = T(n, M, FSi, FS2, RS)A(n, M, FS2, RS), (D.27) 

which follows from (D.18). Differentiating the preceding relation with respect to InM and 
using equation (D.12), we find 

a(M, RS)P(n,M,FSi,RS)T(n,M,FSi,FS2,RS) = dT(n, M FSi, FS2, RS) ^ 

d m M 

+ a(M, RS)T(n, M, FSi, FS2, RS)P(n, M, FS2, RS). (D.28) 

Expanding the preceding equation for e = in powers of a(M, RS), using equation (2.1) 
in the form 

da(/i,RS) ^ _5y-Cfc(RS)a^'+^(;U,RS), co(RS) = 1, ci(RS) = c (D.29) 
dm u ^-^ 

^ k=0 

and after some manipulations, we get the desired formula 

T('=)(n,FSi,FS2),P(°H™)l - A:6TW(n,FSi,FS2) = p('')(n, FSi, RS) - 

k~l 

- P^'^) (n, FS2, RS) + {P^''"'^ {n, FSi , RS)t(') (n, FSi , FS2) - 

- t(') (n, FSi , FS2)P^''"') (n, FS2, RS) + Z6cfc_/(RS)T(') (n, FSi, FS2) } (D.30) 



^^If we use the splitting functions for labeling factorization schemes, then the complete specification of the 
factorization scheme requires also the specification of the corresponding renormalization scheme because the 
relation between the splitting functions and the functions A['^'^-'{x), which define the factorization scheme, 
depends on the renormalization scheme. 
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where k>l. The preceding formula forms a set of equations for T^'^^(n, FSi, FS2), which 
can be solved iteratively. 

The formula for J^5(RSi, RS2, FS) follows from 

a(M, RSi)P(x,M,FS,RSi) = a(M, RS2)P(x, M, J^cS(RSi, RS2, FS), RS2), (D.31) 

which is a consequence of (D.8) and (D.12). Expanding the preceding equation in powers 
of a(M, RS2) (the left hand side is expanded using (B.12)), we obtain 

k 

P(^) {x, -F5(RSi, RS2, FS), RS2) = 4'+/^(RSi, RS2) P^'^ (x, FS, RSi). (D.32) 

1=0 

Contrary to formula (D.30), the preceding relation holds for general e. 
D.5 Transformation of coefficient functions 

The subject of this subsection is the determination of the dependence of coefficient functions 
on the unphysical quantities associated with factorization. An arbitrary structure function 
F(x,Q^) is given as 

F{n, Q2) = c(n, Q^) Bin) (D.33) 

where C(n, Q^) represents the corresponding bare coefficient functions. Using relation 
(D.l), we get 

F(n,Q2) = C(n,Q^M,FS,RS) D(n,M,FS,RS) (D.34) 
where the (finite) coefficient functions C(n, Q^, M, FS, RS) are given by 

C(n,Q^M,FS,RS) = C(n, Q^) B(n, M, FS, RS). (D.35) 

The dependence of coefficient functions on the unphysical quantities associated with the 
factorization procedure then follows from the preceding relation. 
From relation (D.35), we find 

dC(n,Q^M,FS,RS) _ ^ 2^ dB(re, M, FS, RS) _ 
dlnM ~ [n,Q ) dlnM ~ 

= C(n, M, FS, RS) A(n, M, FS, RS)^^^^^^^|^^^^. (D.36) 
Equation (D.12) and A(n, M, FS, RS)B(n, M, FS, RS) = 1 then imply that 
dB(n,M,FS,RS) 



A(n,M, FS,RS) 



dlnM 
dA(n,M,FS,RS) 
dlnM 



B(n, M, FS, RS) = -a(M, RS)P(n, M, FS, RS), (D.37) 



■^^The coefficient functions form a row vector whereas tiie parton distribution functions are represented 
by a column vector, and tiierefore tlie multiplication on the right hand side of (D.33) yields a number (a 
matrix 1x1). 
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which allows us to rewrite formula (D.36) in the form 

Changing the factorization scheme from FSq to FS is described by the following formula 

C(n, Q^, M, FS, RS) = C(n, M, FSq, RS) T(n, M, FSq, FS, RS), (D.39) 

which follows from relations (D.18) and (D.35). 

The coefficient functions C(?i, Q^, M, FS, RS) are (at least in principle) fully calculable 
within the framework of perturbative QCD and therefore can be expanded in powers of 
a(/Li, RS): 

oo 

C(n, Q\ M, FS, RS) = ^ a''{fi, RS) C^''\n, Q\ fi, M, FS, RS) . (D.40) 

A:=0 

Substituting the preceding expansion in formulae (D.38) and (D.39), we obtain the follow- 
ing formulae describing the dependence of the coefficient functions C('^^(n, Q^, /i, M, FS, RS) 
on the factorization scale 



dC('=)(n,Q2 /x,M,FS,RS) . 
a in M = - ^ ^ ^ ' 



1=0 

k-l-l 

E 9tt-LiiM,f^,RS)P^^\n,FS,RS) (D.41) 



X 

m=0 



and scheme 

C(*^)(n, Q^,fi, M, FS, RS) = ^ C(')(n, ^u, M, FSq, RS) x 

1=0 

k-l 

E 5S_„(M,^,RS)T(-)(n,FSo,FS). (D.42) 



X 

m=0 



Relations (D.8) and (D.35) imply that 

C(n,Q^M,FS,RSl) = C(n,Q^M,J^5(RSi,RS2,FS),RS2) . (D.43) 

Hence, the simultaneous change of the factorization scheme and the renormalization scheme 
from {FSo,RSo} to {FS,RS} can be converted to changing the factorization scheme from 
FSq to J^cS(RS, RSo, FS) for the fixed renormalization scheme RSq, similarly as in the case 
of parton distribution functions. 

Formulae (D.41), (D.42) and (B.7) are sufficient for changing all unphysical parameters 
associated with the renormalization and factorization procedure (even in the case if the 
renormalization scheme of the coupling parameter that is employed for expanding the 
coefficient functions is different from the renormalization scheme used for the factorization 
procedure) . 
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D.6 Transformation of hard scattering cross-sections for lepton-hadron colli- 
sions 



According to the parton model, any inclusive cross-section cr{P) (in general differential) 
depending on observables P and describing a lepton-hadron collision is given as^^ 



a{P) = y] [ dxa^{x,P)Di{x). 
Jo 



(D.44) 



The preceding relation expresses the cross-section cr(P) in terms of bare quantities. Using 
relations (D.3), (A. 9), (A. 10) and (A. 11), we get 



(D.45) 



(D.46) 



a(P) = V / dxa,{x,P,M,FS,RS)Di{x,M,FS,RS) 
i Jo 

where the (finite) hard scattering cross-section ai{x, P, M,FS,RS) is given by 
(7i(x,P,M,FS,RS) = V ['dyd,{xy,P)B,,{y,M,FS,RS). 

The inversion of the preceding relation has the form 

ai{x,P) = y] ['dyaj{xy,P,M,FS,RS)A,,{y,M,FS,RS), (D.47) 

which follows from formulae (D.46), (A.9), (A.ll), (A. 12) and (A. 13). 

The dependence of the hard scattering cross-section ai{x, P, M, FS, RS) on the unphys- 
ical quantities associated with the factorization procedure can be investigated in a similar 
way as in the case of the coefficient functions in the preceding subsection. From equations 
(D.46), (D.47) and (A. 12), we find 



dai{x,P,M, FS,RS) 



dlnM 

= E 

jk 



E 



a,{P) 



dBji{M, FS,RS) 
dlnM 



cjfc(P, M, FS, RS) Akj{M, FS, RS) 



dBji{M, FS,RS) 
dlnM 



(x). (D.48) 



Applying equation (D.37) then gives 
dai{x,P,M, FS,RS) 



dlnM 



i(M, RS)V / dyaj{xy,P,M,FS,RS)Pji{y,M,FS,RS). (D.49) 
. ^0 



*The formula (D.33) for structure functions is a particular case of this formula. 
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The formula for changing the factorization scheme from FSq to FS follows from relations 
(D.46), (D.47) and (A.12): 

ai{x,P, M,FS,RS) = 

= Y\aj{P,M,FSo,RS)(^ (^Ajk{M,FSo,RS)*Bki{M,FS,RS)^ (x). (D.50) 

jk 

Taking into account relation (D.18), we can write the preceding equation as 

fTi(x,P,M,FS,RS) = V / dyaj{xy,P,M,FSo,RS)Tji{y,M,FSo,FS,RS). (D.51) 

Expanding the hard scattering cross-section <Tj(a;, P, M, FS, RS) in powers of a{fi, RS): 

oo 

cTi(a;, P, M, FS, RS) = ^a^+''o(/i,RS)cjP(x,P,^,M,FS,RS), ko > (D.52) 

k=0 

and inserting the expansion into equations (D.49) and (D.51), we obtain the formulae 
describing the dependence of the hard scattering cross-sections al''\x, P, fx, M,FS,RS) on 
the factorization scale 



d(7f^(x,P,/i,M, FS,RS) 
dlnM 



E / i E (^2/' ^' ^' RS) X 

k-l-l ^ 

E 5&-Li(^,M,RS)P,(r)(y,FS,RS) (D.53) 



m,=0 

and scheme 



1 r 

af )(x, P, /X, M, FS, RS) = / dy J ^ (xy, P, /x, M, FSq, RS) x 

j •'^ I «=o 

k-l ^ 

X 5] 5K_„(M,^,RS)r,(r)(y,FSo,FS) . (D.54) 

m=0 J 

The preceding formulae (D.53) and (D.54) represent an analogy of formulae (D.41) and 
(D.42) and together with formula (B.7) are sufficient for changing all unphysical quantities 
associated with the renormalization and factorization procedure. 

D.7 Transformation of hard scattering cross-sections for hadron-hadron col- 
lisions 

In the case of a hadron-hadron collision, any inclusive cross-section cr(P) depending on 
observables P can be expressed as 



a{P) = y2 / dxidx2aij{xi,X2,P)Di/H^{xi)Dj/H^{x2). (D.55) 
/o Jo 
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Using relations (D.3), (A. 9), (A. 10) and (A. 11), we can rewrite the preceding formula for 
(t(P) in the form 



f7(P) = V/ / dxid2;2CTij(xi,X2,P,Mi,FSi,RSi,M2,FS2,RS2) X 
Jo Jo 

X A/Hi(^i,Mi,FSi,RSi)Z),/^,(x2,M2,FS2,RS2) (D.56) 

where the hard scattering cross-section cjjj(xi, X2, -P, Mi, FSi, RSi, M2, FS2, RS2) is given 
as 

kl 



o-jj(xi,X2,-P, Mi,FSi,RSi,M2,FS2,RS2) = / / dyidy2aki{xiyi,X2y2, P) x 

kl J 0-^0 

X fifc,(yi,Mi,FSi,RSi)Si,(y2,M2,FS2,RS2). (D.57) 
The preceding relation (D.57) can be written as 

c7i,(xi,X2,P,Mi,FSi,RSi,M2,FS2,RS2) = 

= V/ dy5fc,(xiy,X2,P,M2,FS2,RS2)Sfci(y,Mi,FSi,RSi) (D.58) 
k Jo 

where 

5fc,(xi,X2,P,M2,FS2,RS2) = V / dyaM{xi,X2y,P)Bij{y,M2,FS2,RS2). (D.59) 

/ ^0 

Relation (D.58) has the same form as relation (D.46), and therefore we can immediately 
write the formulae for changing the factorization scale Mi 

dajf (xi, X2, P, Ml, FSi, RSi, M2, FS2, RS2) 



din Ml 
„i ffe-i 



Yl I dy] J^ai5,)(xiy,X2,P,/i,Mi,FSi,RSi,M2,FS2,RS2) X 



, 1=0 

k-l-1 •\ 

E 9tll-liMl,^i,RSl)pi^\y,FSuRSl)\ (D.eo) 



X 

m=0 



and the factorization scheme associated with hadron Hi (from FS\ to FSi) 

^rff (xi , X2 , P, //, Ml , FSi , RSi , M2 , FS2 , RS2) = 
1 ( ^ 

= Y. I dyJ^^7g(xiy,X2,P,/i,Mi,FsS°\RSi,M2,FS2,RS2) 

r Jo y 1=0 



k-l -\ 

X E 5d-™(Mi,/.,RSi)r(r^(y,FSf\FSi) . (D.61) 

m=0 J 
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Similarly, rewriting formula (D.57) in the form 



ai,(xi,X2,P,Mi,FSi,RSi,M2,FS2,RS2) = 

= V/ dyaa{xi,X2y,P,Mi,FSi,RSi)Bij{y,M2,FS2,RS2) (D.62) 

where 

5iz(xi,X2,P,Mi,FSi,RSi) = V / dyakiixiy,X2,P)Bkiiy,Mi,FSi,RSi), (D.63) 

k 

we obtain the formulae describing the change of the factorization scale M2 

dajf (xi, X2, P, /u, Ml, FSi, RSi, M2, FS2, RS2) _ 
dlnM2 ~ 

= -Y. dy<^ J^aJ?(xi,X2y,P,/i,Mi,FSi,RSi,M2,FS2,RS2) X 

k-l-1 ^ 

X 5^ <7&ii_i(M2,/.,RS2)Pi7^(y,FS2,RS2) (D.64) 

m=0 J 

and the factorization scheme associated with hadron H2 (from FS2^^ to FS2) 

^rjf (xi , X2 , P, //, Ml , FSi , RSi , M2 , FS2 , RS2) = 
1 ( ^ 

= Y.I dyJ V4^(xi,X2y,P,/i,Mi,FSi,RSi,M2,FS(°\RS2) X 

r L=0 

k-l ^ 
X E 5S-™(M2,/.,RS2)li7^(y,FSf ,FS2) . (D.65) 

m=0 J 

It is important to point out that formulae (D.60) and (D.61) correspond to the ex- 
pansion in powers of a(;U, RSi) whereas formulae (D.64) and (D.65) correspond to the ex- 
pansion in powers of a(/i, RS2)- As well as in the case of hard scattering cross-sections for 
lepton-hadron collisions, formulae (D.60), (D.61), (D.64), (D.65) and (B.7) are sufficient for 
changing all unphysical parameters associated with the renormalization and factorization 
procedure. 
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